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Abstract 

Experiments on the oscillatory motion of a suspended bar magnet throws light on the damping 
effects acting on the pendulum. The viscous drag offered by air was found the be the main contributor 
for slowing the pendulum down. The nature and magnitude of the damping effects were shown to be 
strongly dependent on the amplitude. 
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I. INTRODUCTION 



The simple pendulum is pedagogically a very important experiment. As a result of it's 
simplicity, it is introduced at high school level. The equation of motion (EOM) of a pendulum 
undergoing simple harmonic motion (SHM) is given as 

d 2 y 9 2 

whose solution is easily derivable and can be taught in a class which has been introduced to 
calculus. The EOM can be modified to account for damping as seen in a real pendulum and 
yet the equation and it's solution remains trivial as^ 

d 2 y fb\dy 2 

~dW + UJ * oW 

y (t) = e' l3t (Acosu't + BsinJt) (1) 

where j3 = (b/2mL) and u' = ^ — f3 2 . However, this approach taken by textbooks over- 
simplifies the complex motion of the pendulum and implies that only the pendulum's amplitude 
attenuates with time. On the contrary along with the amplitude even the oscillation's time 
period varies 2 -, a feature overlooked in classroom physics and carried forward for a long time 
by students. The difficulty in measuring these variations also does not encourage routine 
experimentation in high schools/ undergraduate laboratories. However, with the advent of 
micro-computers such measurements can now be made easily. Most of the experiments reported 
have measured the change in amplitud e 21 ? 14 while examples of measuring variation in time 
period is rare^. Since, both amplitude and the time period varies with successive oscillation, 
one can expect the pendulum's velocity to vary with time at a given position. While Gregory 1 
used knowledge of the oscillation time period to extract information on the pendulum's velocity, 
Avinash Singh et a& used a novel method to estimate the pendulum's velocity. A bar magnet 
was attached to a rigid semicircular aluminum frame of radius 'L' which pivoted about the 
center of the circle such that the bar magnet oscillates through a coil kept at the mean position. 
As the magnet periodically passed through the coil, it generated a series of em/ pulses. The 
arrangement with proper circuitry determined the peak emf. Avinash et a$ approximated the 
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peak em/(f max ) as 




where co> max is the maximum velocity as the bar magnet passed through the mean position. 
This method has it's advantage when one proposes to study the damping effects in a pendu- 
lum. Most of the works studying the variation in oscillation amplitud e 1 18 with time have the 
pendulum's suspension connected to a variable resistance (potentiometer) which introduces 
a sliding friction in the pendulum's motion. Complex mathematics with assumption that all 
damping contributors act independently is then used to filter out information of each con- 
tribution. Wang et al^ used a novel but costly method using Doppler effect to monitor the 
position of the pendulum to study it's damping. Thus, Avinash et al& provides a interesting 
yet cheap method to study the damped pendulum. While they rightly pointed out that several 
parameters of the experiment such as velocity and strength of the magnet and the number of 
turns in the coil can be varied, they did not explicitly discuss them theoretically or study these 
factors experimentally. Hence, in this manuscript, we have furthered the study made in ref 6 
and have tried to address these issues. 



II. EXPERIMENTAL SETUP 



Our pendula was made by suspending a standard bar magnet by a cotton thread. The 
tread was fastened to a small hook drilled into one pole of the bar magnet. The length of the 
bar magnet (21) was 7cm and the cotton thread (L s ) used was 53cm long. A coil of 1000 turns 
was kept near the pendula's mean position at a distance 'd' from the magnet's lower pole (see 
fig 1). The magnetic field at point 'A' is evaluated by 

\x m 



4tt YBA 2 CA 2 \ y ' 

where m is the dipole moment. 'AC and 'BA' can be written in terms of the pendulum's 
position (angle 0) using the cosine law. That is, 

BA 2 = OB 2 + OA 2 - OB.OAcosQ 
AC 2 = OC 2 + OA 2 - OC.OAcosQ 
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where 



O 




FIG. 1: Pendulum with the mass being replaced by a bar magnet. The detecting coil is at 'A'. 
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hence, 



BA 2 = L 2 + (L a + 21 + df - L S .{L S + 21 + d)cosQ 

AC 2 = (L s + 2l) 2 + (L s + 21 + d) 2 - (L s + 21). (L. + 21 + d)cos& 

Based on the assumption that 21 and d are relatively small compared to L s , the higher 
powers of 21 and d can be neglected. Hence, eqn(J2J) can be written as 



B 



47T 



21 



L 2 s (L s + Al + d)(2-cosQ)_ 
The induced emf is proportional to the rate of change in the number of magnetic lines cutting 
the coil. 



dB 
~dt 
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L 2 s (L s + Al + d) 
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Based on this, the respective induced emf can be written as 
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L 2 a (L a + 41 + d) 

where N is the number of turnings in the coil. EqnfjSJ) can be written in a compact form 

sin® 



(3) 
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L 2 s {L s + Al + d) 

Thus, as the distance between the magnet and detecting coil is increased, the induced emf 
decreases. Infact the induced emf is quite weak and is amplified by an op amp circuit. The 
high input impedance of the IC741 opamp ensures that a true measurement of the emf is 
made. To digitalise this analog signal (see fig 1, ref 6) using an Analog to Digital Convertor 
ADC-0809 (see fig 2a), the amplified output is rectified and the peak value is held by charging 
a capacitor. The capacitor is discharged via a large resistance so that it retains the peak value 
till the next peak value arrives. 

We require the ADC to start conversion once the peak value is attained by the capacitor. 
This implies a synchronization between the input em/pulses and the ADC's start of conversion 
(SOC) pulses. To achieve this synchronization it is best to generate the required SOC pulse 
by wave-shaping the input itself. The amplified input after rectification is fed to a comparator 
which compares to +lv. This is to avoid spurious/accidental triggering due to noise. The 
infinite gain results in pulses with sharp edges. The width of these pulses are approximately 
T Q /4 (for our pendula m 390ms). This would be too large for serving as a SOC and hence 
is reduced to a 5/is pulse using a monostable timer made with IC555— . The sequencing 
and synchronization can be understood from the various waveforms shown in fig 2b. The 
designed circuit digitalises the analog emf and on completion sends an EOC to the computer 
or microprocessor kit (in case of a microprocessor this is done through a programmable I/O 
IC8155 chip, details of which can be found in the book by Goankai^) which then reads the 
eight bit data and stores it for retrivial. This project was done using an 8085 microprocessor 
kit. The programme and flowchart used is detailed in the Appendix. 
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FIG. 2: The (a) schematic diagram of the circuit used and the (b) important waveforms at the points 
marked in the circuit. 

The reliability of our circuit can be tested by measuring the maximum emf induced in 
the coil for varying distances 'd'. EqP| shows that the measured maximum emf would be 
directly proportional to £ Q which inturn is inversely proportional to 'd' (see eqn|3j). Fig 3 



6 



4 



3.5 - 



3 



2.5 - 



2 - 



1.5 



1.2 



1.4 



1.6 



1.8 



2 



2.2 



2.4 



2.6 



2.8 



d (in cm) 



FIG. 3: Variation in the maximum induced emf with increasing distance between the coil and the 
magnet. 

shows the variation in the experimentally determined £ max with 'd'. While the inverse nature 
is evident, the value of (L s + Al) as returned by curve fitting eqn© on our data is substantially 
off mark from the actual lengths. This is expected since eqn(@J and (0) are very simplified 
approximations. 

III. VARIATION OF INDUCED EMF WITH INITIAL DISPLACEMENT 
A. While undergoing undamped oscillation 

The velocity of an undamped pendulum undergoing SHM is given as 



where u = y g/ (L s + 21) is the frequency of oscillation and m is the initial displacement given 
to the pendulum. Therefore, the emf induced by our pendulum undergoing undamped SHM 
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would be given as (using eq HJ) 



02 



(2 - cos0) 2 



(6) 



The variation in induced emf with time of an undamped pendulum undergoing SHM is as cal- 
culated using eqn ([S|) is shown in fig (4). The maximum angular displacement used to generate 
the graph using eqnflnj) was 5°. The emf pulse shown in fig 4 is only for half a cycle starting 
from one extreme position to the opposite extreme. As the magnet approaches the coil, the 
flux increases and as it crosses the mean position, the emf is negative since the magnet is 
receding from the coil. Eventhough the velocity (dO/dt) is maximum at the mean position, 
since the variation in flux (d0/dt) is zero, the induced emf is zero as the pendulum passes the 
mean position. 
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FIG. 4: A measure of the induced emf with oscillating angle. The graph was generated using eqn© 
with m = 5°. 



The position of the pendulum when the maximum emf is generated (between < < @ m ) 
can be found as a problem of maxima and minima 
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2 - cosQf [sinQ± (f ) + cosQ (f )] - (2 - cosQ)^Q (f ) 



(2 - cosQf 

(2 - cosQ) [sznQ^ (f ) + cosQ (f )] - ,m 2 9 (f 
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For cases of small angle oscillations eqn(??) reduces to 



dQ ^° 



e ds hr +(1 - e ) hr 



(7) 
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Solving the quadratic equation 



e 4 -(2 + e 2 rt )0 2 + e 2 rt = o 



we have 



6 



6, 



peak 
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(8) 



Since eqn((ZJ) was used to determine position of extrema, the above condition is only valid for 
undamped small angle oscillations. The maxima as per this condition for magnet oscillating 
through G m = 5° occurs at ±0.0617 radians (or ±3.53°). The maximum em/ that is induced, 
hence is (use eqnEJ 

,e m 



Sin- 



V2 



( 2 - cos^) 

Since, these equations and conditions are essentially valid for small angles, 



'^4A el 



(9) 



However, a physical pendulum is prone to damping and hence in the next section we investigate 
as to how the maximum induced emf varies with initial displacement for a damped pendulum. 
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FIG. 5: Variation of induced emf when oscillation is damped (/3 = 0.45s 1 ) is compared with the case 
of no damping (i.e. (3 = 0.0s _1 ). 
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FIG. 6: The variation of induced emf of fig 3 is plotted w.r.t. angular position of the pendulum (i.e. 
0) for the cases (a) (3 = 0.0s- 1 and (b) = 0.45s- 1 . 

B. While undergoing damped oscillation 

We have already stated in our introduction that the damped motion described by eqnl f l j) 
exhibits how the pendulum's oscillation amplitude decreases exponentially with time. The 
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EOM whose solution is given by eqnf^J describes a linear system. The solution can be further 
trivialized without losing any generality as 



6 = ® m e~ pt sin{uj't) 



(10) 



from which the velocity can be calculated as 

dQ 



^ = uu'e m e- pt cos(uu't) - (3e m e~ pt sin(Lj't) 
= Q m e~ pt [uj'cos(uj't) - 0sin(u't)] 



(11) 



Substituting the above expression in eqn(0J) we obtain the relation showing the variation of 
induced em/ with time. This variation is shown in fig 5. It is also clear from the figure that 
the peaks in the induced em/ occurs at ut = (2n + 1)^. Hence, the angles at which maxima 
occur in general is written as 

. e„ 
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(12) 



where tan</> = uj' / (3. Our circuit is designed only to measure peak emfs at n=0,2,4,6 , where 

only the positive solutions of eqn ()12j) would contribute. Using our condition on eqn ljll)) and 
eqn(jlj) we have 
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(13) 



For small angle oscillations eqn (fT3*|) reduces to 



£,peak 
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(14) 



The variation in em/ (seen w.r.t time in fig 5) when viewed w.r.t oscillating angle G shows 
how the peak position decreases (eqlT2"|) as also the amplitude of the maximum induced emf 
decreases (eqn IT4*|) with each half cycle. EqnfjHJ) and eqn (fT4*|) shows that the maximum emf 
induced for damped pendulums undergoing SHM is directly proportional to the square of 
maximum angular displacement given to the pendulum. We have recorded the first maxima 
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FIG. 7: The variation of induced emf with the initial angular displacement of the pendulum. It shows 
the expected parabolic dependence (i.e. Sp- 
reading (i.e. n=0) of the induced emf for various angles upto 40°. The linear relation between 
£ max and 6^ is evident. Before commenting further, it must be recollected that eq(l) is valid 
for small angle oscillations, i.e. for 9 m < 5°, yet a good linearity is obtained till 9 = 40°. 

Experimental data for O m > 45° deviate markedly from this linear trend. Remember 
eqn (jl2j) was obtained with the assumption that the pendula's motion is described by eqn (jl()|) . 
This equation describes the motion of a pendulum oscillating in a viscous medium with small 
velocity. It would be shown in the next section that the pendulum's velocity is quite apprecia- 
ble for 9 m > 45° and hence it's motion is not described as in eqn lfTIH) . explaining the departure 
for linearity. 

IV. RESULTS AND OBSERVATIONS 

Our prelimary measurements are in good correspondence with commonly known notions 
and hence we proceed to investigate further the nature of damping in our pendula. It should 
be noted that the amount of damping and it's nature are strongly pendula dependent and all 
results reported here are specific to our experiment and can not be taken as general. We have 
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FIG. 8: The variation in peak induced emf measured with each oscillation is shown for initial dis- 
placements (e m ) (a) 5°, (b) 30°, (c) 55° and (d) 65°. 



recorded the maxima in induced emfs for 80 oscillations for various initial displacements. Since 
for each oscillations, we get two positive maxima in induced emf figure 8 shows the variation 
in maxima reading of induced emf for 160 peaks. 

A general expression ae~ bn was fitted to the data of £ max w.r.t n using a standard and freely 
available curve fitting software called "Curxpt v3.0". Good fits were obtained for oscillations 
set by initial displacements upto 40°. The exponential fall in emf is indicative of the rate of 
loss of energy from the oscillating system. It indicates the loss to follow the relation 

dE 

— — oc —E 
dt 

or 

dt 

This indicates that the velocity is low and hence the damping/resistive force acting on the 
pendulum is proportional to the velocity^. Figure 9 shows the variation of the decay constant 
'b' (of ae~ bn ) with respect to the maximum displacement (0 m ) given to the pendulum. The 
graph indicates that as G m increases, the velocity with which the pendulum moves increases 
with which the damping constant increases. 
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FIG. 9: The variation decay constant (b of general equation ae bn ) of the exponentially decaying 
region. The continuous line is the parabolic fit for the data points (0. 00056 v / ©m)- 



FIG. 10: The variation decay in 'b', a measure of decay in the early oscillations when pendulum was 
set in motion with displacements > 45°. 



For displacement angles beyond 45° (0 m > 45°), eventhough visual examination of the 
curves in fig 8 suggests an exponential damping, the data points do not fit to an exponential 
fall relation. A detailed examination suggests a more complex process is taking place with 
initial damping being sharper. Infact the initial 25-30 data points fit to a/n b . The data 
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beyond this fit to the exponential fall equation. The a/n b fit corresponds to the damping force 
being proportional to higher powers of velocity (v 7 , where 7 > 1) and in turn the rate of energy 
loss also being proportional to higher terms of energy. That is, the rate of energy loss for our 
pendulum set into oscillations with a displacement angle > 45° is given as 

dE i+b 
— — = —aE b 
dt 

Figure 10 plots 'b' versus G m . The power term (b) is being treated as a measure of damping 
and is consistent with the results of fig 9, i.e. as initial displacement increases the damping 
becomes large with a proportionality to the pendulum's velocity. 

The resistive force being proportional to higher powers of velcoity has been reported earlier 
a l S0 ^ddM^l^ilL, A system is reported to have a constant friction (7 = 0) or a linear dependence 
of velocity (7 = 1) or a quadratic depenence of velocity (7 = 2). Corresponding to which 
the pendulum's amplitude decays linearly, exponentially and inverse power decay respectively 
with time. It hence maybe concluded that for our pendulum sent into motion by initial 
displacements > 45°, the damping force is proportional to v 7 where 7 > 1. 



V. CONCLUSION 



A simple experiment of setting a suspended bar magnet into oscillations, is a rich source of 
information. Not only does it give exposure to Faraday's induction law and a basic understand- 
ing of induced emf 's dependence on angle of oscillation, it enables us to study the damping 
effects on the pendulum. This method is better than previously used methods since the mea- 
suring technique does not introduce additional contributions to damping. When the oscillation 
imparted to the pendulum is very large, the damping effect is also strong with the damping 
force being proportional to v 7 , where 7 > 1 and V is the pendulum's velocity. This brings 
down the oscillation amplitude of the pendulum and it's velocity. As the velocity becomes 
low, the resistive force acting on the pendulum changes it's nature and becomes proportional 
to V. Considering the rich information obtained from the experiment and the simplicity of 
the experiment, it allows the method to be easily implemented as a routine experiment in 
undergraduate laboratories. 
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Appendix 



Table 1. Program used to collect data. 
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FIG. 11: Flowchart. 
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